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SUMMARY 

An existing computer code that determines the flow path for an axial-flow compressor either for a 
given number of stages or for a given overall pressure ratio was modified for use in air-breathing engine 
conceptual design studies. This code uses a rapid approximate design methodology that is based on isen- 
tropic simple radial equilibrium. Calculations are performed at constant-span-fraction locations from tip 
to hub. Energy addition per stage is controlled by specifying the maximum allowable values for several 
aerodynamic design parameters. 

New modeling was introduced to the code to overcome perceived limitations. Specific changes 
included variable rather than constant tip radius, flow path inclination added to the continuity equation, 
input of mass flow rate directly rather than indirectly as inlet axial velocity, solution for the exact value 
of overall pressure ratio rather than for any value that met or exceeded it, and internal computation of 
efficiency rather than the use of input values. The modified code was shown to be capable of computing 
efficiencies that are compatible with those of five multistage compressors and one fan that were tested 
experimentally. 

This report serves as a users manual for the revised code, which is named CSPAN, an acronym for 
Compressor SPanline AN alysis. The modeling modifications, including two internal loss correlations, are 
presented. Program input and output are described. A sample case for a multistage compressor is 
included. 


INTRODUCTION 

Performing engine studies requires the capability to produce conceptual designs of the components 
in order to determine geometry, performance, and weight. One major component of air-breathing turbine 
engines is the compressor. The typical compressor “design” code enables a study of the interrelationship 
of the number of stages, the flow path radii, the gas velocities, the flow angles, and the resultant varia- 
tion of compressor efficiency. A computer code capable of performing this function in a rapid approximate 
manner was selected as being consistent with the needs of engine conceptual design. This code (ref. 1), 
which is based on isentropic simple radial equilibrium, was one of four compressor analysis codes devel- 
oped under NASA contract about 25 years ago. The other three codes, two design (refs. 2 and 3) and one 
off-design (ref. 4), were streamline analyses accounting for full radial equilibrium. 

An evaluation of the reference 1 code indicated several limitations to its usefulness for the con- 
ceptual sizing of engine compressors. The design was restricted to a constant tip radius, and flowpath 
inclination was not accounted for in the continuity calculation. In addition, the user was unable to 
directly specify a mass flow rate or an exact overall pressure ratio. Finally, stage efficiency had to be 
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estimated beforehand because it was a required input, and there was no provision for inlet guide vane 
loss. Consequently, the code was modified to overcome these deficiencies. 


This report serves as a users manual for the revised code, which is named CSPAN, an acronym for 
Compressor SPanline ANalysis. The modeling changes are presented herein. Computed efficiencies are 
compared with test results from five multistage compressors and one fan. Program input and output are 
described. A sample case for a multistage compressor is included, 

SYMBOLS 

B coefficient in tangential velocity equation (eq.(2)), (ft)(in.)/sec 
b axial distance between calculation stations, ft 
C coefficient in tangential velocity equation (eq.(2)), ft/sec 

D coefficient in tangential velocity equation (eq.(2)), ft/ (sec) (in.) 

D diffusion factor 

E coefficient in tangential velocity equation (eq.(2)), ft /(sec) (in. 2 ) 
f fraction of span height 

g gravitational constant, 32.17 (lbm)(ft)/(lbf)(sec 2 ) 

H enthalpy, Btu/lb 

J conversion constant, 778 (ft)(lb)/Btu 

k loss coefficient multiplier 

M Mach number 

n stage number 

P pressure, psi 

R radius, ft 

T temperature, °R 

U blade speed, ft/sec 

V velocity, ft/sec 

w mass flow rate, lb/sec 

z loss function 

a angle of inclination of flow or endwalls in meridional plane, deg 
/9 flow angle on conical blade-to-blade surface, deg 
rj efficiency 

p density, lb/ft 3 

a solidity 

(j pressure-loss coefficient 
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Subscripts: 


avg 

average 

cor 

corrected 

e 

equivalent 

ex 

exit 

g 

geometric 

H 

hub 

id 

ideal 

in 

inlet 

inp 

input 

ns 

normal shock 

ov 

overall 

P 

poly tropic 

pr 

profile 

R 

rotor 

S 

stator 

s 

static 

sh 

shock 

stg 

stage 

T 

tip 

t 

total 

z 

meridional component 

9 

tangential component 

1 

rotor inlet 

2 

rotor exit 

2D 

two dimensional 

3 

stator exit 

3D 

three dimensional 

Superscripts: 

t 

relative 

* 

uncorrected value 



METHOD OF ANALYSIS 


The computer code of reference 1 was developed for making parametric studies of advanced multi- 
stage axial-flow compressors. This code determines the meridional flow path for given design specifica- 
tions of either the number of stages or the overall pressure ratio. Such a flow path is illustrated in 
figure 1(a), where the stage station locations and some of the geometry variables are defined. A typical 
stage velocity diagram including the symbols for all velocities and angles is shown in figure 1(b). 

The flow- physics model and the solution procedure are described in detail in reference 1. These are 
briefly summarized herein, and then the revisions made to improve the code’s usefulness for air-breathing 
engine conceptual design studies are described. 


Flow-Physics Model 


The radial equation of motion is based on isentropic simple radial equilibrium (no effects of stream- 
line slope and curvature and radially constant entropy). Consequently, at each axial station 


dH t V, d(RV,) 

gj - — + 

dR R dR 


dV. 

dR 


( 1 ) 


where the variation of tangential velocity with radius is 
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Equations (l) and (2) along with continuity 


w = 2 jt f RT p V Z R dR 
jr h 


( 2 ) 


( 3 ) 


and stage energy addition 


gj AH t = A(UV,) 


( 4 ) 


provide the basic flow-physics model. 

Stage energy addition is determined either by specifying the tangential velocities (eq. (2)) at both 
the rotor inlet and exit or by specifying the tangential velocity at the rotor inlet, the axial velocity ratio 
across the rotor tip, and a maximum value for the rotor tip diffusion factor, 
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R,T 




2a V, 


( 5 ) 
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from which rotor-exit tangential velocity is obtained. In either case, the energy addition is reduced if 
limit values are exceeded for stator-inlet hub Mach number, stator hub diffusion factor, or rotor-exit hub 
relative flow angle. Energy addition is then related to rotor and stage pressure ratios by input values of 
rotor and stage polytropic efficiencies, respectively. 

Computation proceeds stage by stage until either a given number of stages is reached or a given 
overall pressure ratio is met or exceeded. The flow path geometry evolves from the given inlet tip radius 
and continuity and the specified limits for hub and tip ramp angles. Tip radius normally remains con- 
stant and hub radius is computed. If the hub ramp angle limit is exceeded, the tip radius is reduced so 
that the hub ramp angle is at its limit value. 


Code Revisions 

The code of reference 1 was evaluated for application to the conceptual design of compressors for 
air-breathing engines. These conceptual designs serve as the basis for estimating engine weight and per- 
formance. It was found that the code’s usefulness could be improved by revising some of its physical 
modeling and solution procedures. The basic flow-physics modeling and energy-addition modeling were 
retained as previously described except for the inclusion of flow path inclination in the continuity equa- 
tion (eq. (3)) as discussed later in this section. All modeling changes are described in this section. 

Tip radius . — With the original methodology the compressor design was executed with a constant 
tip radius unless a hub ramp angle constraint was exceeded. Many advanced designs require tip radius 
reductions in order to provide adequate blade height at the exit. Therefore, an input was added to the 
code that allowed direct specification of the tip radius change across each blade row. The hub ramp 
angle limit was kept as a constraint in order to avoid excessive wall slope and if exceeded would result in 
further tip radius reduction. 

Continuity . — In the original model the velocity is expressed in terms of two components 

V 2 = V 2 + V 2 (6) 

z 9 

The V z term is referred to as “axial” velocity and is used as such for continuity (eq. (3)) even though 
there can be significant inclination in parts of the flow path (especially in the hub region of the inlet 
stages). As used in equation (6), V z is actually a meridional velocity (the resultant of axial and radial 
components). Considering V z as a meridional velocity herein, the continuity equation is revised to 

w = 2?r f R ' r pV z cos a R dR (7) 

jr h 


The streamline angle of inclination a is estimated as the slope of the constant-span-fraction line through 
the previous blade row. 

Mass flow rate . — The original code required as input a value for inlet tip axial velocity, which was 
then used to compute a mass flow rate for the compressor. Because mass flow rate is usually specified for 
an engine, it is preferable to use it rather than inlet axial velocity as the input. Therefore, the calculation 
procedure was modified to enable a mass flow rate to be input and an iteration to be performed to find 
the value of inlet tip axial velocity that provided the given mass flow rate. 
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Overall pressure ratio . — With the original code the analysis proceeded stage by stage either until a 
given number of stages was reached or until a given pressure ratio was met or exceeded. This procedure 
did not provide an exact solution for the given overall pressure ratio. Therefore, the code was modified 
so that convergence to the exact specified pressure ratio can be achieved. This is done by maintaining 
the number of stages constant once the specified pressure ratio is exceeded and then reducing the max- 
imum allowable rotor-blade loading (i.e., the rotor-exit tangential velocity) for all stages until the desired 
overall pressure ratio is obtained. 

Rotor hub turning . — One of the aerodynamic constraints for a design is the amount of turning at 
the rotor hub. In the original code this was specified by an input value for rotor-exit hub relative flow 
angle. An appropriate limit value for this variable, however, is a function of the specific design. In order 
to generalize this constraint, it was changed to a direct input specification of rotor hub turning angle. 

Inlet guide vane loss . — -Computations are begun at the rotor inlet and therefore do not include the 
inlet guide vanes. Although a tangential velocity distribution can be specified at the first-rotor inlet (to 
simulate the exit flow from the inlet guide vanes), there is no way to include the inlet guide vane pressure 
loss as part of the compressor overall pressure ratio. In order to avoid accounting for this loss by arti- 
- ficially reducing compressor inlet pressure and increasing overall pressure ratio, a pressure-loss fraction 
AP/P for the inlet guide vanes was added to the input and included as part of the overall pressure ratio. 

Internal loss correlations .— In the original code, rotor and stage polytropic efficiencies, which were 
assumed to be radially constant, had to be specified as input for each stage. As a consequence there was 
no assurance of compatibility of the input values with the stage loading or the blade-element aerodynam- 
ics. Therefore, two loss correlations were added to the code for optional use: one for stage polytropic 
efficiency and the other for blade-element, pressure-loss coefficient. For engine conceptual design studies 
the stage polytropic efficiency correlation is recommended. 

Stage poly tropic efficiency: Axial-flow compressor efficiency correlations (unpublished) that are 
used at NASA Lewis Research Center for engine cycle studies are presented in figure 2. Stage polytropic 
efficiency is plotted against stage pressure ratio in figure 2(a) for current- and advanced- technology axial- 
flow compressor stages. The advanced-technology curve represents the improvement that is anticipated 
over the next 10 to 20 years. Shown in figure 2(b) is the efficiency correction for small compressors (i.e., 
for low values of corrected mass flow rate). The curves of figure 2 were fit with the following equations: 
For current-technology stages the equation for pressure ratios of 2 or less is 




0.054795 



0.25337 



+ 1.1477 


( 8 ) 


A linear extrapolation is used beyond a stage pressure ratio of 2 to yield 


^P.atg 


-0.03419 



+ 0.9285 


( 9 ) 
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For advanced-technology stages the equation for pressure ratios of 2 or less is 


Vstg 


0.047322 



0.21668 



1.1241 


( 10 ) 


An extrapolation beyond a stage pressure ratio of 2 yields 




-0.027392 



+ 0.93478 


( 11 ) 


For corrected mass flow rates of less than 10 lb/sec, efficiency is reduced to account for size effects 
(clearances, surface finish, etc.). The size correction for corrected mass flow rates between 1.5 and 
10 lb/sec is 


Aij p gtg = 0.56826xl0~ 3 w cor - 0.62224xl0 -3 


0.050603 


w 


cor 


( 12 ) 


and linear extrapolation for flow rates below 1.5 lb/sec yields 


A»7 P)8tg = 0.01767 w cor - 0.06 


(13) 


The corrected flow rate in equations (12) and (13) is defined as 

/ T in /518.7 


w 


W„ 


(14) 




After correcting for size, the resultant value of stage polytropic efficiency can be adjusted by an input loss 
multiplier as follows: 


Vstg - 1 - k i 


inp 


1 (Vatg + 


(15) 


Assuming that two-thirds of the stage loss occurs in the rotor, the rotor polytropic efficiency is 


^P)R 1 g ( 1 ’Wig) 


( 16 ) 


This approximation, which is based on experience, affects only the hub radius at rotor exits. Equa- 
tions (8) to (16) were incorporated into the code as one optional method for computing efficiency. 
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Blade-element pressure loss: Another method for estimating efficiency is through the use of loss 
coefficients that are based on blade-element aerodynamics. Rotor and stator loss coefficients are defined 
as 


W R 




Pt,2 “ Pt,3 
Pt,2 ~ Pg,2 


(17a) 


(17b) 


Each blade-element loss coefficient will be composed of profile loss and shock loss components. 

The profile loss coefficients are obtained from two-dimensional cascade data with a correction for 
three-dimensional effects. Two-dimensional loss coefficient data as a function of diffusion factor (eq. (5)) 
were obtained from reference 5 and extrapolated by using trends from reference 2. These loss data were 
then fit with the curve 


Z 2D “ 


W 2D cos £ex 
2a 


= 0.0065 + 0.0050566 D + 0.027721 D 4 ' 7773 


(18) 


The data, extrapolation, and curve fit are presented in figure 3. The default curve-fit coefficients used for 
equation (18) can be replaced through program input by alternative values that match other sources of 
loss data. 

The loss function z was then modified for three-dimensional effects on the basis of the loss data 
used in reference 6 as follows: 


Z 3D ~ ^3D Z 2D 


(19) 


The three-dimensional correction 


k 3 D is based on the fraction of span height from the tip 


Ro> - R 
f = 

R'p ~ Rjj 


For the outer 30 percent of span (f < 0.3) 


k 3D - 1.87 - 2.9f 


and for the inner 30 percent of span (f > 0.7) 


k 3D = 2f - 0.4 


( 20 ) 


(21a) 


(21b) 
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There is no correction for the center 40 percent of blade span. 


The profile loss function was increased by 50 percent on the basis of experimental compressor per- 
formance (see next subsection) and can be arbitrarily adjusted further by an input loss multiplier. 


V = 15k mp z 3D 


( 22 ) 


The profile loss coefficient is then 


2az : 


oj = 
P r 


pr 


cos ft 


(23) 


The shock loss coefficient is based on the methodology of reference 7, wherein the passage-shock 
loss is taken as the normal-shock loss from the arithmetic average of the inlet and the estimated peak 
suction-surface Mach numbers. The authors of reference 6 state that a shock loss so determined is 
excessive, and they recommend that the normal-shock loss be divided by the square of the average Mach 
number. Therefore, the rotor shock loss coefficient is 


(ij = 

ns 


i - p. /p. , 

t,ns' t,l 


(24) 


1 - P 




and 






sh 


/ \2 
M 

avg/ 


(25) 


Although a shock also can theoretically occur in the stator, aerodynamic constraints that are imposed on 
the design normally prevent this. 


The blade-element overall loss coefficient is then 

; pr T w sh 


GJ ~ (jJ + 0J 

w Qy w nr w « 


(26) 


and the blade-row exit total pressure is determined from this loss coefficient. 

Comparison with experimental performance: In order to test the loss models, calculated perfor- 
mance was compared with experimental performance for five multistage compressors (refs. 8 to 12) and 
one fan (ref. 13). The overall design features of these machines are presented in table I along with esti- 
mates of the design-point efficiencies that are based on test data. 

The data reported in references 8 to 13 are for research or early-development compressors and 
therefore are not representative of the design-point efficiency levels that are achievable in developed 
engines. Consequently, for the purposes of this comparison the measured design-point efficiencies of the 
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referenced compressors were adjusted to a projected design-point value. The projected value was taken as 
the maximum efficiency either at design speed or, if there appeared to be a significant mismatch at design 
speed, at a somewhat lower speed. This adjustment procedure is largely subjective, but there appear to 
be no published data on the performance of fully developed current- and advanced-technology compressors. 

The referenced compressors were modeled for the CSPAN program, which was run using both loss 
models. The design-point efficiency comparison of calculated and measured efficiencies is presented in fig- 
ure 4. For the poly tropic efficiency correlation the calculated and measured efficiencies were within one 
point of each other for all five multistage compressors and were within two points for the fan. For the 
pressure-loss coefficient correlation the multiplier of 1.5 introduced into equation (22) was based on the 
best overall comparison between calculated and measured values. With this value, four of the compres- 
sors compared to within about one point, but the calculated efficiency of the fan was more than three 
points low. These loss models seem to yield reasonable estimates for multistage compressors, but some 
adjustment in the input loss multiplier k inp may be needed for fans. 

DESCRIPTION OF INPUT AND OUTPUT 

This section presents a detailed description of input and output for program CSPAN. Included 
with the input and output is a sample case for a five-stage transonic compressor. 


Input 


The input, which is read on unit 05, consists of a title line and one NAMELIST dataset. Input for 
the sample case is presented in table II. The title, which is printed as a heading on the output file, can 
contain up to 71 characters located anywhere in columns 2 through 72 on the title line. A title, even if it 
is left blank, must be the first record of the input data. 


The physical data and option switches are input in data sets having the NAMELIST name NAME. The 
variables that compose NAME are defined herein along with units and default values. They are pre- 
sented in order as general inputs, inlet inputs, rotor inputs, and stator inputs. 


General: 

CP 

MW 

GAM 

RCLIM 

NSLIM 

N 

ICV 


specific heat of working fluid, Btu/(lb)(°R) 
molecular weight of working fluid, lb/ (lb mol) 
specific heat ratio 

limit value for overall pressure ratio 

limit value for number of stages 

number of calculation locations from tip to hub 

pressure ratio convergence switch (default — 1) 

0— accepts overall pressure ratio equal to or greater than RCLIM as a 
solution 

1 — converges to overall pressure ratio equal to RCLIM 
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IPR1 debug output switch (default = 0) 

0 — no debug output 

1 — minimum debug output 

2 — extensive debug output 

IPR2 station output switch (default = 0) 

0 — output printed for all radial locations 

1 — output printed for tip and hub only 

DPPIGV inlet guide vane total-pressure loss fraction (default = 0.0) 

WK loss multiplier, equation (15) or (22) (default = 1.0) 

IIT technology-level indicator for polytropic efficiency (default = 1) 

1 — current technology, equations (8) and (9) 

2 — advanced technology, equations (10) and (11) 

AZ0 constant term in pressure-loss coefficient correlation, equation (18) (default = 0.0065) 

AZ1 coefficient of linear term in pressure-loss coefficient correlation, equation (18) 

(default = 0.0050566) 

AZ2 coefficient of exponential term in pressure-loss coefficient correlation, equation (18) 

(default = 0.027721) 

AZ3 exponent in pressure-loss coefficient correlation, equation (18) (default = 4.7773) 

Inlet: 

TTI inlet total temperature, °R 

PTI inlet total pressure, psi 

RTIP1I tip radius at first-rotor inlet, in. 

UTIP1I blade speed at first-rotor inlet, ft /sec 

RHORTl hub/tip radius ratio at first-rotor inlet 

VZTIPO inlet axial velocity or mass flow rate specifier 

>0 — VZTIPO is the axial velocity at first-rotor inlet, ft /sec 
<0 — [VZTIPO I is the mass flow rate, Ib/sec 

DTIP1 tip blockage factor at first-rotor inlet (default = 1.0) 

DH1 hub blockage factor at first-rotor inlet (default = 1.0) 

B1 coefficient B for equation (2) at first-rotor inlet, (ft)(in.)/sec 

Cl coefficient C for equation (2) at first-rotor inlet, ft/sec (default = 0.0) 

D1 coefficient D for equation (2) at first-rotor inlet, ft/(sec)(in.) (default = 0.0) 

El coefficient E for equation (2) at first-rotor inlet, ft /(sec) (in. 2 ) (default = 0.0) 

Rotor: Each variable requires NSLEM values. 

RT20Tl(I) ratio of exit tip radius to inlet tip radius for each rotor (default = 1.0) 

VT20Tl(I) ratio of exit tip meridional velocity to inlet tip meridional velocity for each rotor 
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NPRI(I) efficiency specifier for each rotor 

> 0.0— input value is rotor polytropic efficiency 
_ o.O — polytropic efficiency correlation is used 

= —1.0 — pressure-loss coefficient correlation is used 

SRTIP(I) rotor tip solidity 

ARO(I) rotor aspect ratio (based on axial chord) 

DTIP2(I) tip blockage factor at rotor exit (default = 1.0) 

DH2(I) hub blockage factor at rotor exit (default = 1.0) 

ARHD(I) rotor hub ramp angle limit, deg 

ARTD(I) rotor tip ramp angle limit, deg 

DRT(I) rotor tip diffusion factor maximum value 

BO(I) rotor-exit tangential velocity specifier 

= 0.0— rotor tip exit tangential velocity determined from rotor tip 
diffusion factor by equation (5) 

> 0.0 — coefficient B for equation (2) at rotor exit, (ft) (in.) /sec 

C2(I) coefficient C for equation (2) at rotor exit, ft/sec (default — 0.0) 

D2(I) coefficient D for equation (2) at rotor exit, ft/(sec)(in.) (default — 0.0) 

E2(I) coefficient E for equation (2) at rotor exit, ft/(sec)(in. 2 ) (default = 0.0) 

BPSD(I) limit value for rotor hub turning, deg 

Stator: Each variable requires NSLEM values. 

RT30T2(I) ratio of exit tip radius to inlet tip radius for each stator (default = 1.0) 
VT30T2(I) ratio of exit tip meridional velocity to inlet tip meridional velocity for each stator 

NPSI(I) efficiency specifier for each stage 

> 0.0 — input value is stage polytropic efficiency 
= 0.0 — polytropic efficiency correlation is used 

= —1.0 — pressure-loss coefficient correlation is used 

SSH(I) stator hub solidity 

ASO(I) stator aspect ratio (based on axial chord) 

DTIP3(I) tip blockage factor at stator exit (default = 1.0) 

DH3(I) hub blockage factor at stator exit (default = 1.0) 

ASHD(I) stator hub ramp angle limit, deg 

ASTD(I) stator tip ramp angle limit, deg 

DSH(I) limit value for stator hub diffusion factor 

MSH(I) limit value for stator hub inlet Mach number 

B3(I) coefficient B for equation (2) at stator exit, (ft)(in.)/sec 

C3(I) coefficient C for equation (2) at stator exit, ft/sec (default = 0.0) 
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D3(I) 

E3(I) 


coefficient D for equation (2) at stator exit, ft/(sec)(in.) (default = 0.0) 
coefficient E for equation (2) at stator exit, ft/(sec)(in. 2 ) (default — 0.0) 


Output 

Program output consists of a main output file written to unit 06 and, if applicable, a brief 
pressure-ratio convergence file written to unit 08. The main output presents either the results of a 
successful design calculation or an error message indicating the nature of the failure to find a solution 
that is consistent with the design specifications. 

Outputs corresponding to the sample input of table II are presented in tables III and IV. The 
pressure-ratio convergence output, shown in table III, is most useful when sent to the terminal so that a 
convergence problem can be immediately detected and the computation halted. Convergence problems, 
however, have not occurred for any of the six cases tested (table I). As shown in table III, convergence to 
a pressure ratio of 5 required four iterations. Shown in the output are the number of stages followed by 
one line for each iteration displaying the rotor tip diffusion reduction factor (DRTK), the compressor 
pressure ratio (CPR), and the compressor adiabatic efficiency (EFF). 

The main output is presented in table IV. For brevity in displaying the output, calculations were 
performed at only three radial locations (N = 3), and only the data for stages 1 and 5 are included along 
with the overall and inlet information. The first line of output in table IV is the title; it is followed by 
identification of the loss model used for this case. Then, the general inputs and the inlet inputs are 
printed. The values displayed are clearly identified. 

The next output line in table IV states that one of the aerodynamic limits, in this case the stator 
hub Mach number, was exceeded in the next stage, which in this case is the first stage. As a result, the 
rotor tip diffusion factor was reduced from its maximum allowable value until the Mach number limit 
was just satisfied. The consequence of this is a reduction in stage pressure ratio. 

The next block of output in table IV is the data for stage 1. This includes the rotor input, the 
stator input, the stage performance and geometry, and the detailed aerodynamic results at the rotor inlet, 
the rotor exit, and the stator exit. Note that the rotor and stator input sections include the rotor and 
stage polytropic efficiencies, respectively, which were^Jetermined in this case from the internal correla- 
tion. Under stage output data are the overall values of pressure ratio, temperature ratio, and adiabatic 
efficiency followed by the stage values, which are the same since this is the first stage. Also displayed are 
the rotor and stator tip and hub radii, the axial lengths, and the tip and hub ramp angles. Finally, at 
each of the three axial stations for the stage are presented among other parameters, the temperatures and 
pressures, the absolute and relative velocities, the absolute and relative flow angles, the diffusion factors, 
and the loss coefficients at each of the radial calculation locations. 

The stage data format is identical for each stage; therefore, the “STAGE DATA” output for 
stages 2, 3, and 4 were omitted from table IV. Shown next is the data for stage 5. The “STAGE 
OUTPUT DATA” show that the overall pressure ratio of 5 has been achieved with an overall efficiency 
of 0.8775. For this constant-tip-radius (10 in.) design, the hub radius increased from 5.0 in. at the first- 
rotor inlet to 8.5 in. at the last-stator exit. The last line of output states that the specified overall pres- 
sure ratio has been achieved. If the specified pressure ratio had not been achieved, the last-line message 
would be that the maximum number of stages had been reached. 
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SUMMARY OF RESULTS 


An existing computer code that determines the flow path for an axial-flow compressor either for a 
given number of stages or for a given overall pressure ratio was selected for use in air-breathing engine 
conceptual design studies. This code uses a rapid approximate design methodology that is based on isen- 
tropic simple radial equilibrium. Calculations are performed at a number of constant-span-fraction loca- 
tions from tip to hub at each blade-row inlet and exit. Energy addition per stage is controlled by a 
maximum allowable value for the rotor tip diffusion factor, which is reduced if limit values are exceeded 
for the stator hub inlet Mach number, the stator hub diffusion factor, or the rotor hub turning angle. 

This code was modified to make it easier to use for the conceptual study of engine compressors. 

The rapid approximate design methodology was retained and new modeling was introduced to overcome 
perceived limitations. The limitations and associated modifications were as follows: 

1. Unless the given hub ramp angle limit was exceeded, the tip radius had remained constant. A 
tip radius change for each blade row can now be specified through input. 

2. The throughflow component of velocity had been assumed to be purely axial for the continuity 
calculation. An internally computed flow path inclination angle was added to the continuity equation. 

3. Mass flow rate had been computed from an input value of inlet axial velocity. A new algorithm 
allows mass flow rate to be input and the proper value of inlet axial velocity to be found by iteration. 

4. Any case wherein the overall pressure ratio met or exceeded the specified value had been taken 
to be an acceptable solution. Convergence to the exact value of overall pressure ratio can now be 
achieved. 

5. Rotor and stage polytropic efficiencies had to be input for each stage. Two internal loss cor- 
relations were added to the code as alternative ways to specify performance: One is for the stage poly- 
tropic efficiency and the other is for the blade-element, pressure-loss coefficient. For engine conceptual 
design studies the stage polytropic efficiency correlation is recommended. 

6. There had been no provision to include inlet guide vane pressure loss as part of the overall pres- 
sure ratio. An input pressure-loss fraction for the inlet guide vanes has been added. 

The modified code was tested by comparison with five multistage compressors and one fan for 
which experimental data were available. The computed performance was found to be compatible with the 
test results. 

This report serves as a users manual for the modified code, which is named CSPAN, an acronym 
for Compressor SPanline AN alysis. Program input and output are described. A sample case for a multi- 
stage compressor is included. 
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TABLE I.— COMPRESSORS USED FOR LOSS MODEL EVALUATION 


Reference 

StageB 

Overall 

pressure 

ratio 

Tip 

radius, 

in. 

Tip 
speed, 
ft /sec 

Measured 

efficiency 

8 

8 

10.3 

10.0 


0.87 

9 

5 

5.0 

10.0 


.87 

10 

3 

4.5 

10,1 

1 

.86 

U 

10 

23.0 

13.8 

1 

.86 

12 

10 

14.0 

11.7 

1 

.86 

13 

1 

1.65 

41.4 

K9i 

.89 


TABLE II.— SAMPLE INPUT 


NACA 5 STAGE TRANSONIC COMPRESSOR 
&NAME 

CP= . 24 ,MW=29 . , GAM=1 . 4 , RCLIM^S . 0 , NSLIM=5 , N=3 , 

TTI=518 . 7,PTI=14.7,RTIP1I=10. ,UTIP1I=1100 . ,RH0RT1= . 5 , VZTIPO=-67 . 5, B1=0 . 0 , 
VT20T1= . 918 , .885, .916, . 909 , .916, NPRI=5*Q . 0 ,SRTIP=0 . 98 , 1 . 17 , 1 . 30, 1 . 1* , 0 . 99 , 
ARO-2 . ,1.52,1 .11 ,0 .99,0 .92, DRT=5*. 45,50-5*0. 0,BPSB=5*45, , 

VT30T2=1 .094,1 .071,1 . 053 , 1 . 057 , 0 . 992 ,NPSI=5*0 . 0 ,SSH=1 .8 , 1 . 9 ,1 . 6 , 1 . 5 , 1 . 4 , 
A$0=2, 15,1 .63,1 .24,1.01 , 0 . 88 , DSH-5* . 55 , MSH=5* . 75 , B3=5#0 . 0 , 

SEND 


TABLE III.— CONVERGENCE OUTPUT FOR SAMPLE CASE 


STAGES- 5 

DRTK- 1.00000000 
DRTK- 0.916593611 
DRTK= 0.924198091 
DRTK= 0.92444610 6 


CPR- 5.45497990 
CPR= 4.95435619 
CPR= 4.99855804 
CPR 3 5.00001240 


EFF= 0.873050213 
EFF= 0.877941251 
EFF= 0.877498388 
EFF= 0.877484918 
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TABLE IV— MAIN OUTPUT FOR SAMPLE CASE 


NACA 5 STAGE TRANSONIC COMPRESSOR 


LOSS MODEL: INTERNAL CORRELATION FOR STAGE POLYTROPIC EFFICIENCY 


*** INLET INPUT DATA XX* 


NO. RAD. 
STATIONS 

3 

NUMBER 

STAGES 

5 

SP. HEAT 
(BTU/(LB-R)) 

0.2A00 

MOL. WT. 
(MOLES) 

29.0000 

RATIO OF 
SP. HEAT 

1 .4000 

IN. TOT. TEMP. 
(DEG. R) 

518.7000 

IN. TOT. PR. 
(PSI) 

1A.7000 

MASS A V6. TOT. 
PR. RATIO 

5.0000 

IG 

DEL 

0.0 



XX * ROT 

OR INLET 

INPUT D 

A T A xxx 





TIP RADIUS 
(INCHES) 

10.0000 

TIP WHEEL SPEED 
(FT/SEC) 

1100 . 0000 

HUB TO TIP 
RADIUS RATIO 

0.5000 

MASS FLOW 
(LB'/SEC) 

67.5Q0Q 

TIP BLOCKAGE 
FACTOR 

1.0000 

HUB BLOCKAGE 
FACTOR 

1 .0000 





COEFFICIENTS IN TANGENTIAL VELOCITY 

EQUATION 






B 

C 

D 

E 






0 . 0000 

0.0000 

0.0000 

0 .0000 





XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXKXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 


XXX STATOR HUB MACH NO. LIMIT VIOLATED X** 
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TABLE IV. — Continued. 


******** STAGE DATA ******** 
STAGE NO. 1 


*** ROTOR INPUT DATA *** 


merid vel. 

RATIO 

POLYTROPIC 

EFFICIENCY 

SOLIDITY 
AT TIP 

ASPECT 

RATIO 

TIP 

BLOCKAGE 

FACTOR 

HUB 

BLOCKAGE 

FACTOR 

MAX ANGLE 
HUB TAPER 
(DEGREES) 

MAX ANGLE 
TIP TAPER 
(DEGREES) 


0.9180 

0.9308 

0.9800 

2.0000 

1.0000 

1.0000 

40.000 

-20.000 


MAX ROTOR 
DIF. FACTOR 

MAX. TURNING 
ANGLE ROTOR HUB 
(DEGREES) 

TIP 

RADIUS 

RATIO 


COEFFICIENTS IN TANGENTIAL VELOCITY 
BCD 

EQUATION 

E 


0.4160 

45.0000 

1. 000 


0.0000 

0.0000 

0,0000 

0 . 0000 





*** STATOR 

INPUT DATA *** 





MERID 

VELOCITY 

RATIO 

STAGE 

POLYTROPIC 

EFFICIENCY 

SOLIDITY 
AT HUB 

ASPECT 

RATIO 

TIP 

BLOCKAGE 

FACTOR 

HUB 

BLOCKAGE 

FACTOR 

MAX ANGLE 
HUB TAPER 
(DEGREES) 

MAX ANGLE 
TIP TAPER 
(DEGREES) 


1.0940 

0.8963 

1.8000 

2.1500 

1.0000 

1.0000 

40.0000 

-20.0000 


MAX. STATOR 
DIF. FACTOR 

MAX HUB INLET 
MACH NUMBER 

TIP 

RADIUS 

RATIO 


COEFFICIENTS IN TANGENTIAL VELOCITY 
BCD 

EQUATION 

E 


0.5500 

0.7500 

1.000 


0.0000 

0.0000 

0.0000 

0.0000 




***--**--*** 

STAGE 0 

UTPUT DATA 

***--**--*** 







MASS FLOW (LB/SEC) * 67.501 





OVERALL 
MASS AVE. 
PR. RATIO 

OVERALL 
MASS AVE. 
TEMP. RATIO 

OVERALL 
MASS AVE. 
EFFICIENCY 

HASS AVE. 
PRESSURE 
RATIO 

MASS AVE. 
TEMPERATURE 
RATIO 

MASS AVE. 
EFFICIENCY 

ROTOR 

ASPECT 

RATIO 

STATOR 

ASPECT 

RATIO 

OVER 

MASS 

POLY 

1.4423 

1.1238 

0.8907 

1.4423 

1.1238 

0.8907 

2.0000 

2.1500 

0. 

ROTOR TIP 
RAD. i"G 
(INCHES) 

ROTOR HUB 
RAD. 1-G 
(INCHES) 

ROTOR TIP 
RAD. 2-G 
(INCHES) 

ROTOR HUB 
RAD. 2-G 
(INCHES) 

STATOR TIP 
RAD. 3-G 
(INCHES) 

STATOR HUB 
RAD. 3-G 
(INCHES) 

ROTOR PROJ. 
LENGTH 
(INCHES) 

STATOR PROJ. 
LENGTH 
(INCHES) 


10.0000 

5.0000 

10.0000 

5.8071 

10.0000 

6.4608 

2.5000 

1.9502 




ROTOR TIP 
RAMP ANGLE 
(DEGREES) 

ROTOR HUB 
RAMP ANGLE 
(DEGREES) 

STATOR TIP 
RAMP ANGLE 
(DEGREES) 

STATOR HUB 
RAMP ANGLE 
(DEGREES) 






0.0000 

17.8924 

0.0000 

18.5305 








***--**--*** R 

0 T 0 R 

I N L E 

T OUTPUT DATA *** 

--**--*** 





STA 

NO. 

RADIUS 

-E 

(IN) 

WHEEL 

SPEED 

(FT/SEC) 

MERID 

VEL. 

(FT/SEC) 

TANGENT . 
VEL. 

(FT/SEC) 

ABS. 

VEL. 

(FT/SEC) 

REL. 

VEL. 

(FT/SEC) 

ABS. REL. TOTAL 

AIR ANG. AIR ANG. TEMP. 
(DEG) (DEG) (DEG R) 

TOTAL 

PRESS. 

(PSI) 

REL. 

MACH 

NO. 

ABS. 

MACH 

NO. 

SHOCK 

LOSS 

COEFF 

TOTAL 

LOSS 

COEFF 

ROTO 

DIF. 

FACT 

1 

2 

3 

10.000 

7.500 

5.000 

1100.000 

825.000 

550.000 

637.702 

637.702 

637.702 

0.000 

0.000 

0.000 

637.702 

637.702 

637.702 

1271.480 

1042.731 

842.118 

0.000 59.898 518.700 
0.000 52.297 518.700 
0.000 40.777 518.700 

14.700 

14.700 

14.700 

1.179 

0.967 

0.781 

0.591 

0.591 

0.591 

0 . 000 
0.000 
0.000 

0.048 

0.065 

0.090 

0.393 

0.456 

0.472 




#**-- 

■**--*** 

ROTOR 

: E x i i 

' OUTPUT DATA ***- 

-**--*** 





STA 

NO. 

RADIUS 

-E 

(IN) 

WHEEL 

SPEED 

(FT/SEC) 

MERID 

VEL. 

(FT/SEC) 

TANGENT. 

VEL. 

(FT/SEC) 

ABS. 

VEL. 

(FT/SEC) 

REL . 
VEL. 

(FT/SEC) 

ABS. REL. TOTAL 

AIR ANG. AIR ANG. TEMP. 
(DEG) (DEG) (DEG R) 

TOTAL 

PRESS. 

(PSI) 

SPANL 

INCL 

(DEG) 

REL. 

MACH 

NO. 

ABS. 

MACH 

NO. 

LOSS 

FUNC 

STA! 

REACT 

1 

2 

3 

10.000 

7.904 

5.807 

1100.000 

869.391 

638.702 

585.410 

585.410 

585.410 

350.745 

443.781 

603.991 

682.442 

734.606 

841.136 

950.835 

723.774 

586.443 

30.928 51.999 582.931 
37.165 36.018 582.931 
45.895 3.401 582.931 

21.503 

21.503 

21.503 

0.000 

9.170 

17.892 

0.832 

0.637 

0.523 

0.597 

0.646 

0.750 

0.015 

0.021 

0.025 

0.92: 

0.82i 

0.61! 




***--**--*** 5 

; T A T 0 

R E X I 

T OUTPUT DATA ***--**--*** 





STA 

NO. 

RADIUS 

-E 

(IN) 


MERID 

VEL. 

(FT/SEC) 

TANGENT, 

VEL. 

(FT/SEC) 

ABS. 

VEL. 

(FT/SEC) 

SPANL. 

INCL. 

(DEG) 

ABS. TOTAL SHOCK 

AIR ANG. PRESS LOSS 

(DEG) (PSI) COEFF 

TOTAL 

LOSS 

COEFF 

STATOR 

DIF. 

FACTOR 

MERID 

MACH 

NO. 

ABS . 
MACH 
NO. 

LOSS 

FUNC 


1 

2 

3 

10.000 

8.230 

6.461 


640.438 

640.439 
640.439 

0.000 

0.000 

0.000 

640.438 

640.439 
640.439 

0.000 

9.514 

18.530 

0.000 21.201 0.000 

0.000 21.201 0.000 

0.000 21.201 0.000 

0.065 

0.057 

0.045 

0.294 

0.349 

0.438 

0.558 

0.558 

0,558 

0.558 

0.558 

0.558 

0.030 

0.021 

0.013 
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TABLE IV.— Concluded. 


XXXXXXXX STAGE DATA XXXXXXXX 
STAGE HO. 5 

xxx ROTOR INPUT DATA XXX 


HER ID VEL. 
RATIO 

POLYTROPIC 

EFFICIENCY 

SOLIDITY 
AT TIP 

ASPECT 

RATIO 

TIP 

BLOCKAGE 

FACTOR 

HUB 

BLOCKAGE 

FACTOR 

MAX ANGLE 
HUB TAPER 
(DEGREES) 

MAX ANGLE 
TIP TAPER 
(DEGREES) 


0.9160 

0.9416 

0.9900 

0.9200 

1.0000 

1.0000 

40 .000 

-20. 0C0 


MAX ROTOR 
DIF. FACTOR 

HAX . TURNING 
ANGLE ROTOR HUB 
(DEGREES) 

TIP 

RADIUS 

RATIO 


COEFFICIENTS IN TANGENTIAL VELOCITY 
BCD 

EQUATION 

E 


0.4160 

45.0000 

1.000 


0 . 0000 

0.0000 

0.0000 

0 . OOOO 





XXX STATOR 

INPUT DATA XXX 





MERID 

VELOCITY 

RATIO 

STAGE 

POLYTROPIC 

EFFICIENCY 

SOLIDITY 
AT HUB 

ASPECT 

RATIO 

TIP 

BLOCKAGE 

FACTOR 

HUB 

BLOCKAGE 

FACTOR 

MAX ANGLE 
HUB TAPER 
(DEGREES) 

MAX ANGLE 
TIP TAPER 
(DEGREES) 


0.9920 

0.9123 

1.4000 

0.8800 

X .0000 

1.0000 

40 .0000 

-20.0000 




TIP 


COEFFICIENTS IN TANGENTIAL VELOCITY 

EQUATION 


MAX. STATOR 
DIF. FACTOR 

MAX HUB INLET 
MACH NUMBER 

RADIUS 

RATIO 


B 

C 

D 

E 


0.5500 

0.7500 

1 .000 


a . oooo 

0 .0000 

0 .0000 

0. oooo 




XXX- -XX- -XXX 

STAGE Q 

UTPUT DATA 

xxx - -xx- - xxx 







MASS FLOW (LB/SEC) - 67.501 





OVERALL 
HASS AVE. 
PR. RATIO 

OVERALL 
HASS AVE. 
TEMP. RATIO 

OVERALL 
MASS AVE. 
EFFICIENCY 

HASS AVE- 
PRESSURE 
RATIO 

MASS AVE. 
TEMPERATURE 
RATIO 

MASS AVE. 
EFFICIENCY 

ROTOR 

ASPECT 

RATIO 

STATOR 

ASPECT 

RATIO 

OVERALL 
MASS AVE 
POLY EFF. 

5.0000 

1 .6653 

0.8775 

1.2874 

1.0823 

0.9091 

0.9200 

0.8800 

0.9016 

ROTOR TIP 
RAD. 1-G 
( INCHES! 

ROTOR hub 
RAD. 1-G 
(INCHES) 

ROTOR TIP 
RAD. 2-0 
(INCHES) 

ROTOR HUB 
RAD. 2-6 
(INCHES) 

STATOR TIP 
RAD. 3-6 
(INCHES) 

STATOR HUB 
RAD. 3-G 
(INCHES) 

ROTOR PROJ. 
LENGTH 
(INCHES) 

STATOR PROJ. 
LENGTH 
(INCHES) 


10.0000 

8.3072 

10.0000 

8.4424 

10.0000 

8.4819 

1.8400 

1.7700 




ROTOR TIP 
RAMP ANGLE 
(DEGREES) 

ROTOR HUB 
RAMP ANGIE 
(DEGREES) 

STATOR TIP 
RAMP ANGLE 
(DEGREES) 

STATOR HUB 
RAMP ANGLE 
(DEGREES) 






0.0000 

4.2029 

0.0000 

1.2772 





xxx- -XX- -xxx ROTOR INLET OUTPUT DATA XXX--XX--XXX 


• 

ST A 
NO. 

RADIUS 

-E 

(IN) 

WHEEL 
SPEED 
(FT /SEC) 

MERID 

VEL. 

(FT/SEC) 

TANGENT. 

VEL. 

(FT/SEC) 

ABS. 

VEL. 

(FT/SEC) 

REL . 
VEL . 

(FT/SEC) 

ABS. 

AIR ANG. 
(DEG) 

REL. TOTAL 

AIR ANG. TEMP. 
(DEG) (DEG R) 

TOTAL 

PRESS. 

(PSI) 

REL. 

MACH 

NO. 

ABS. 
MACH 
NO . 

SHOCK 

LOSS 

COEFF 

TOTAL 

LOSS 

COEFF 

ROTOR 

DIF. 

FACTOR 


1 

Z 

3 

10.000 

9.154 

8.307 

1100.000 

1006.897 

913.795 

562.565 

562.562 

562.562 

0 .000 
0.000 
0.000 

562.565 

562.562 

562.562 

1235.506 
1153.394 
1073 . 077 

0.000 

0.000 

0.000 

62.914 79B.105 
60.807 798.105 
58.382 798.105 

57.093 

57.093 

57.093 

0.908 

0.847 

0.7B8 

0.413 

0.413 

0.413 

0.000 
0 . 000 
0.000 

0.039 

0.043 

0.048 

0.416 

0.451 

0.490 





xxx- - 

XX- -xxx 

ROTOR 

! EXIT 

D U T F 

» U T DATA xxx- 

-xx -xxx 






STA 

NO. 

RADIUS 

-E 

(IN) 

WHEEL 

SPEED 

(FT/SEC) 

MERID 

VEL. 

(FT/SEC) 

TANGENT. 

VEL. 

(FT/SEC) 

ABS . 
VEL . 

(FT/SEC) 

REl . 
VEL. 

(FT/SEC) 

ABS. 

AIR ANG . 
(DEG) 

REL. TOTAL 

AIR ANG. TEMP. 
(DEG) (DEG R) 

TOTAL 

PRESS. 

(PSD 

SPANL 

INCL 

(DEG) 

REL . 
MACH 
NO. 

ABS . 
MACH 
NO. 

LOSS 

FUNC 

STAGE 

REACTION 


1 

2 

3 

10.000 

9.221 

8.442 

1100. 000 
1014.334 
928.668 

515.309 

515.307 

515.307 

358.787 

389.088 

424.980 

627.910 

645.701 

667.944 

902.740 

810,230 

720.585 

34.848 

37.055 

39.513 

55.192 863.809 
58.506 863.809 
44.347 863.809 

74.097 

74.097 

74.097 

0.000 

2.104 

4.203 

0 . 639 
0.574 
0.512 

0.445 

0.458 

0.474 

0.011 
0 . 013 
0.015 

0.901 
0.873 
0 836 





XXX- - XX - - xxx 

S T A T 0 

R E X I 

T OUT 

PUT DATA xxx- -xx - -xxx 






STA 

NO. 

RADIUS 

-E 

(IN) 


MERID 
VEL . 

( FT /SEC ) 

TANGENT . 
VEL . 

(FT/SEC) 

ABS. 
VEL . 

(FT/SEC) 

SPANL . 
INCL . 
(DEG) 

ABS. 

AIR ANG. 
(DEG) 

TOTAL SHOCK 

PRESS LOSS 

(PSD COEFF 

TOTAL 

LOSS 

COEFF 

STATOR 

DIF. 

FACTOR 

MERID 

MACH 

NO. 

ABS. 

MACH 

NO. 

LOSS 

FUNC 


c 

fi 

1 

2 

3 

10.000 

9.241 

8.482 


511.187 

511.184 

511.184 

0.000 

0.000 

0.000 

511.187 
511 . 184 
511.184 

0.000 

0.639 

1.277 

0 . ODD 
0.000 
0.000 

73.500 0.000 
73.500 0.000 
73.500 0.000 

0.064 

0.060 

0.056 

0 . 427 
0 . 443 
0.462 

0.360 

0.360 

0.360 

0.360 

0.360 

0.360 

0.027 
0. 023 
0.020 






xxx OVERALL PRESSURE RATIO LIMIT 

HAS BEEN 

REACHED -- GO TO NEW DATA 

xxx 
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Stag© 1 
exit 


Stage 2 
exit 


Stage 3 
exit 


Last -stage 
exit 



v z,1 = Vi — ' U, V e ,3 

(b) Stage velocity diagram. 

Figure 1 .—Schematic presentation of symbols. 


20 




Technology 




(a) Effect of stage pressure ratio. 

Figure 2. — Axial -flow compressor efficiency. 


(b) Size correction. 



Figure 3. — Variation of loss function with diffusion factor. 



Measured efficiency 


Figure 4. — Comparison of measured and calculated 
efficiencies. 
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